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INTRODUCTION 


This  semiannual  report  Includes  the  research  results  achieved 
to  date  at  S-Cubed  (S^)  under  Contract  No.  F49620-81-C-0093.  The 
main  objectives  of  this  contract  are: 

1.  To  provide  support  to  Mission  Research  Corporation  In 
their  work  on  seismic  excitation  of  atmospheric  acoustic 
waves,  and 

2.  To  assess  the  effects  of  block  motions  near  burled 
explosions  on  the  seismic  radiation  field  excited  by  such 
events,  by  numerical  simulations  and  by  analysis. 

This  report  Is  organized  In  three  sections.  The  first  one 
provides  an  executive  summary  of  our  results  to  date,  and  outlines 
the  conclusions  reached  so  far  by  S-Cubed  based  on  these  results. 

The  second  one  focuses  on  the  block  motion  problem,  for  which 
specific  objectives  are  listed  below: 

Task  81:  Formulate  boundary  conditions  for  slipping  joints. 

Task  B2:  Perform  two-dimensional  axi symmetric  calculations  to 
compare  with  experimental  results,  as  well  as 
predictions  of  the  linear  theory,  and  to  assess  the 
results  in  preparation  for  three-dimensional 
simulations. 

Task  B3:  To  perform  three-dimensional  simulations  incorporating 
free  surface  interaction,  gravity,  joint  orientation 
and  impedance  contrast. 

Task  B4:  To  prepare  synthetic  surface  wave  radiation  patterns 
from  such  simulations. 

Task  B5:  To  compare  numerical  and  analytical  models  and 
assess  the  limitations  of  the  latter. 

Task  B6:  To  compare  the  results  with  those  of  a  conventional 
tectonic  release  model  and  assess  their  effects  on 
yield  estimates. 
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The  last  section  pertains  to  research  on  the  coupling  between 
seismic  and  acoustic  waves.  The  corresponding  objectives  are: 

Task  Al:  To  assess  the  uncertainties  In  the  methodology 

chosen  to  evaluate  such  coupling. 

Task  A2:  To  review  the  literature  applicable  to  shallow 

earthquakes  of  magnitude  M  -  4. 

Task  A3:  To  gather  near  field  data  for  such  an  event. 

Task  A4:  To  provide  Mission  Research  Corporation  with 

selsmologlcal  data  for  simulation  of  Ionospheric 
disturbances  due  to  small  earthquakes  and  nuclear 
explosions. 

This  report  represents  the  current  status  of  our 
Investigations  of  these  two  problems,  and  focuses  on  the  objectives 
which  have  been  reached  during  the  first  contract  period. 
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I.  EXECUTIVE  SUMMARY 
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This  report  covers  progress  accomplished  at  S-Cubed  (S  )  on 
research  under  Contract  Number  F49620-81-C-0093. 

The  goals  of  this  research  fall  into  two  main  categories;  (1) 
modeling  the  effects  of  block  motions  in  the  immediate  vicinity  of 
an  underground  nuclear  explosion  on  the  elastodynamic  seismic 
radiation  associated  with  such  events,  and  (2)  evaluation  and 
modeling  of  the  coupling  of  seismic  radiation  to  acoustic  waves  in 
the  atmosphere  and  the  ionosphere. 

Progress  on  these  two  aspects  of  the  research  for  the  period 
September  1,  1981  to  January  31,  1982  can  be  summarized  as  follows. 

1.1  BLOCK  MOTION  MODELING 

We  have  conducted  a  review  of  the  literature  on  boundary 
conditions  at  joints  in  a  geologic  material  (Task  B-l).  Such 
boundary  conditions  fall  into  three  categories:  (1)  linear  boundary 
conditions  which  are  amenable  to  analytic  treatment  of  the  problem, 
and  correspond,  as  a  class,  to  viscous  behavior  of  the  medium;  (2) 
nonlinear  boundary  conditions  which  are  rate-independent,  of  which 
Coulomb  friction  is  the  most  obvious  example;  and  (3)  nonlinear 

boundary  conditions  of  the  type  studied  by  Dieterich  (1979),  which 
include  an  explicit  strain  rate  dependence  and  which  thus  introduce 
an  implicit  time  scale  into  the  problem. 

The  weight  of  laboratory  and  field  evidence  favors  nonlinear 
behavior  for  realistic  models  of  joints  near  an  explosion.  This,  in 
turn,  means  that  simulations  must  be  conducted  numerically.  Linear 
boundary  conditions  should  not  be  Ignored,  however,  in  view  of  the 
fact  that  analytical  Investigations  permit  a  more  general 

understanding  of  the  Interaction  of  joints  with  the  outgoing 

radiation  field  than  a  necessarily  limited  set  of  numerical 

calculations.  Consequently,  our  attack  proceeded  along  the 
following  lines: 
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We  have  conducted  two-dimensional  numerical  simulations  (Task 
8-2)  using  an  axisymmetric  geometry,  for  both  viscous  and  Coulomb 
boundary  conditions.  The  results  have  been  evaluated  against  field 
observations  associated  with  the  MIGHTY  EPIC  test,  as  summarized  by 
Bache,  et_al_.  (1979).  From  these  simulations  we  conclude  that  block 
motions  can  be  modeled  realistically  using  a  simple  Coulomb 
criterion.  Block  motion  is  found  to  occur  in  the  immediate  vicinity 
of  the  shot  point  (inside  the  elastic  radius)  and  to  decrease 
rapidly  with  distance.  These  results  are  both  in  qualitative  and 
quantitative  agreement  with  observations.  Viscous  boundary 
conditions  yield  estimates  of  block  motions  which  are  in  qualitative 
agreement  with  the  estimates  obtained  for  nonlinear  3.C.  However, 
viscous  behavior  predicts  slippage,  and  thus  perturbations  to 
outgoing  waves,  even  for  joints  which  are  far  removed  from  the  shot 
point,  a  feature  which  disagrees  with  observations. 

We  have  compared  the  numerical  results  with  analytical  models 
in  the  case  of  a  viscous  joint  of  unlimited  extent  (Task  B-5).  For 
this  purpose  we  have  generalized  and  corrected  the  analytical 
treatment  of  Salvado  and  Minster  (1980).  We  are  now  in  a  position 
to  define  an  equivalent  source  of  seismic  radiation  which  includes 
the  effects  of  block  motions.  For  linear  behavior,  which  can  be 
modeled  analytically,  this  permits  verification  of  the  numerical 
results,  and  therefore  validation  of  the  numerical  models.  Such 
comparisons  have  been  highly  favorable  and  give  us  confidence  that 
the  approach  is  a  valid  one. 

In  addition,  given  the  results  of  an  arbitrary  numerical 
simulation,  which  might  Include  nonlinear  behavior,  the  analytical 
treatment  also  yields  an  equivalent  source  in  that  case.  This 
greatly  simplifies  the  problem  of  calculating  far-field  seismic 
waves  in  the  case  of  axisymmetric  geometry  (Task  B-4). 
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During  this  contract  period,  we  have  confined  the  study  of  the 
far-field  radiation  to  analysis  of  the  effective  source  spectrum, 
for  cylindrically  symmetric  geometries.  The  equivalent  source  is 
represented  in  terms  of  a  multipolar  expansion,  and  thus  lends 
itself  to  very  convenient  comparison  with  other  causes  of  anomalous 
radiation  (e.g.,  tectonic  release  (Task  &-6)).  Our  results  to-date 
indicate  that  a  passive  mechanism  such  as  block  motion  is  very 
inefficient  in  terms  of  its  ability  to  modify  long  period 
radiation.  In  particular,  we  conclude  that  for  the  axi symmetric 
situation,  block  motion  cannot  explain  the  reversal  of  teleseismic 
20  second  Rayleigh  waves  or  the  generation  of  long  period  Love 
waves.  On  the  other  hand,  spectral  perturbations  at  body  wave 
periods  can  be  significant.  We  are  in  the  process  of  computing 
synthetic  seismograms  (Task  B-4)  in  order  to  quantify  these  effects 
and  gauge  their  observability. 

The  first  portion  of  this  report  provides  a  more  detailed 
description  and  discussion  of  our  results  to-date. 

1.2  ACOUSTIC  COUPLING 

The  main  object  of  this  aspect  of  the  research  is  to  provide 
Mission  Research  Company  (MRC)  with  seismo logical  data  and 
analytical  expertise  to  simulate  ionospheric  disturbances  caused  by 
small  earthquakes  and  underground  nuclear  explosions.  Tasks  A-l  and 
A-2  pertained  to  our  assessment  of  uncertainties  in  the  methodology 
used  for  coupling  seismic  and  acoustic  waves,  and  a  review  of  the 
literature  applicable  to  shallow  earthquakes  of  about  magnitude 
4.0.  These  tasks  have  been  completed  and  the  results  are  summarized 
below  in  Section  III. 

We  have  also  collected  near-field  ground  motion  recordings  of 
several  events  of  magnitude  ~  4.0  (Task  A-3).  For  this  purpose  we 
have  chosen  seven  well-recorded  aftershocks  of  the  Oroville 
aftershock  sequence  and  obtained  near-field  records  from  the 


5 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


U.S.  Geological  Survey  (USGS).  Furthermore,  Me  also  have  obtained 
fault  plane  solutions  and  fault  parameter  estimates  for  these  events 
(J.  Boatwright,  USGS,  personal  communication).  We  are  currently  in 
the  process  of  verifying  that  a  simple  source  model  embedded  in  an 
elastic  half-space  can  adequately  account  for  the  features  of 
near-field  ground  motion  relevant  to  acoustic  coupling.  When  this 
Is  accomplished,  we  shall  be  in  a  position  to  provide  Mission 
Research  Company  with  a  cost-  and  time-effective  means  of  estimating 
ground  motion  from  such  events  at  a  large  number  of  points  at  the 
earth's  surface. 
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II.  EFFECTS  OF  BLOCK  MOTIONS  ON  THE  RADIATION  FIELD 
OF  BURIED  EXPLOSIONS 


2.1  INTRODUCTION 

The  problem  of  anomalous  radiation  from  a  burled  explosion  Is 
ultimately  connected  to  the  question  of  whether  tectonic  stress  Is 
present  or  not.  The  release  of  tectonic  stress  may  take  place  in 
the  form  of  an  actual  earthquake  triggered  by  the  explosion  (e.g., 
Akl  and  Tsai,  1972),  or  In  the  form  of  prestress  relaxation  due  to 
the  creation  of  the  cavity  (e.g.,  Archambeau,  1972;  1973).  In 
either  case,  the  anomalous  component  has  the  radiation  pattern  of  a 
quadrupole  or  a  double  couple,  and  analytical  approximate  solutions 
are  available,  both  in  the  far  field  and  in  the  near  field  (e.g., 
Archambeau,  1968,  1972;  Minster,  1973,  1979;  Minster  and  Suteau, 
1977;  Barker  and  Minster,  1980).  Exact  solutions  are  also  available 
and  have  recently  been  reviewed  by  Stevens  (1980). 

Other  sources  of  anomalous  radiation  also  operate  in  the 
vicinity  of  the  shot  point.  Of  particular  interest  Is  the  class  of 
phenomena  collectively  described  as  "block  motion"  (Bache  and 
Lambert,  1976;  Pratt  and  Rawson,  1977).  It  was  shown  in  previous 
work  (Salvado  and  Minster,  1980)  that  block  motion  due  to 
imperfectly  bonded  interfaces  may  contribute  to  the  generation  of 
anomalous  radiation.  This  mechanism  does  not  require  the  presence 
of  tectonic  prestress.  It  generalizes  the  mode  conversion  mechanism 
proposed  by  Akl  and  Larner  (1970),  and  operates  in  a  fashion  similar 
to  the  cracking  mechanism  of  Kiss  linger,  et  aK  (1961).  We  have 
shown  that  this  mechanism  is  reasonably  efficient  only  if  certain 
constraints  are  met,  but  there  is  reason  to  believe  that  conditions 
may  be  sufficiently  favorable  near  the  shot  point. 

The  connection  between  near  shot  point  block  motion  and 
selsmological  observations  has  been  discussed  and  reviewed  by  Bache 
and  Lambert  (1976).  Specific  estimates  of  block  motions  from 
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sel sinological  observations  Mere  derived  by  Bache,  et  aj_.  (1979)  in 
the  case  of  the  MIGHTY  EPIC  and  DIABLO  HAWK  events.  These  and 
similar  investigations  tend  to  adopt  the  point  of  vieM  that  block 
motion  is  most  conveniently  treated  by  introducing  an  additional, 
spurious  source  field  to  the  primary  source  field.  In  the  present 
study,  we  adopt  the  point  of  view  that  the  primary  source  field  is 
propagated  through  a  set  of  imperfectly  bonded  interfaces  and  loses 
energy  to  alternate  modes  of  wave  motion.  Block  motion  is  then  only 
an  epiphenomenon  and  Its  analysis  is  a  consequence  of  solving  the 
wave  propagation  problem. 

One  key  ingredient  in  the  analysis  of  block  motion  phenomena 
is  the  set  of  boundary  conditions  Imposed  along  the  joint.  We 
provide  below  a  brief  review  of  laboratory  and  field  evidence 
bearing  on  this  problem,  followed  by  the  description  of  a  set  of 
axisymmetric  numerical  calculations  showing  the  effects  of  several 
choices  of  boundary  conditions.  We  then  use  an  analytical  approach 
to  define  an  equivalent  seismic  source  and  obtain  preliminary 
estimates  of  the  effects  of  block  motion  on  teleseismic  radiation 
fields. 


2.2  BOUNDARY  CONDITIONS  AT  A  JOINT  IN  A  ROCK  MASS 

We  can  classify  possible  boundary  conditions  for  rock  joints 
as  (1)  linear,  (2)  nonlinear  and  rate-independent,  or  (3)  nonlinear 
and  rate-dependent.  Upon  review  of  the  relevant  literature,  we 
conclude  the  following: 

1.  The  linear  boundary  condition  is  an  inappropriate  model 
for  frictional  sliding  on  rock  joints. 

2.  A  rate- Independent  Coulomb-type  model,  with  slip 
hardening  to  some  peak  strength  followed  by  slip 
weakening  to  some  kinetic  friction  level,  Is  compatible 
with  the  gross  features  of  laboratory  observations. 


8 


systems.  science  and  software 


3.  Rate-dependent  effects  have  been  well  documented  In  the 
laboratory,  but  only  for  relatively  low  slip  rates. 
Evidence  suggests  that  the  simpler  rate-independent  model 
Is  adequate  at  the  very  high  slip  rates  of  Interest  to  us 
In  the  block  motion  problem. 

Our  preferred  friction  model  Is  summarized  in  Figure  1.  It  Is 
based  largely  on  observations  reviewed  by  Barton  (1976).  The  model 
is  rate-independent,  and  the  residual  shear  stress  is  simply 
proportional  to  the  effective  normal  stress,  afl,  with  coefficient 
of  kinetic  friction  equal  to  0.77.  The  latter  coefficient  has 
been  found  to  be  quite  Insensitive  to  rock  type  (Byerlee,  1968). 
The  value  of  peak  strength  is  a  fairly  strong  function  of  normal 
stress,  as  described  by  Barton. 

2.3  TW0-DIf€NSI0NAL  MODELS 
2.3.1  Introduction 

In  this  section,  we  present  numerical  results  obtained  from 
two-dimensional  ( ax 1 symmetric)  simulations  of  block  motion.  These 
simulations  were  performed  using  a  finite  difference  method.  They 
incorporated  a  nonlinear  friction  law  for  rock  joints  based  on 
laboratory  measurements  as  reviewed  in  the  last  section. 

This  part  of  the  study  has  two  main  objectives.  First,  we 
wish  to  validate  the  accuracy  of  the  numerical  method  for  treating 
displacement  discontinuities  across  joints.  This  was  done  by 
applying  the  numerical  method  to  a  linear  joint  model  and  comparing 
the  numerical  solution  with  the  analytic  solution  of  Salvado  and 
Minster  (1980).  The  second  objective  is  to  establish  the 
plausibility  of  modeling  block  motion  by  the  physical  mechanism  of 
dry  friction  using  frictional  properties  determined  from  laboratory 
samples.  This  was  done  by  comparing  a  nonlinear  simulation  with 
field  observations  of  block  motion  for  the  MIGHTY  EPIC  event. 


9 


systems,  science  and  so rrwAme 


2 


an  tan  [20  log-jg  (200/an)  +  37.5°],  an  <  200  MPa 
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2  *  on  tan  [37.5°] 

Joint  constitutive  model  used  in  nonlinear  simulation  of  block 
moti on . 
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2.3.2  Linear  Simulation 

The  friction  law  considered  here  is  that  treated  analytically 
by  Salvado  and  Minster  (1980).  This  solution  is  examined  in  greater 
detail  In  Section  2.4.3.  It  Is  assumed  that  shear  traction  on  the 
crack  plane,  _r_,  is  proportional  to  the  particle  velocity 
discontinuity,  s: 


In  Equation  (1),  u  is  the  shear  modulus,  a  is  the  shear  wave  speed, 
and  4  is  a  dimensionless  parameter  which  lies  in  the  range  0  to  1, 
and  which  determines  the  strength  of  the  bonding  across  the  joint. 
We  will  treat  the  problem  of  an  explosive  point  source  near  a  rock 
discontinuity  governed  by  Equation  (1).  Note  that  in  the  static 
limit,  the  solution  to  this  linear  problem  will  be  independent  of  4, 
for  4>  <  1.  In  the  numerical  calculations,  we  will  treat  the  case 
*  *  0,  1.  e.,  a  perfectly  lubricated  joint. 

The  problem  geometry  is  shown  in  Figure  2.  The  distance  from 
the  source  point  to  the  fault  is  h,  and  r  is  radial  distance  from 
the  symmetry  axis  (a  line  normal  to  the  fault  through  the  source 
point).  The  source  strength  can  be  described  by  its  reduced 
displacement  potential  4*  (t) ,  with  static  valued,*,.  The  source  time 
history  was  assumed  to  have  the  form  of  a  smoothed  step,  i.e., 

f(t')  -  foe  (.5  -  .5  cos  *  t'/0.6) 

where  t*  is  the  dimensionless  time  t8/h.  The  computed  slip  time 
history  at  three  different  ranges  is  shown  in  Figure  3,  where  It  is 
compared  with  the  analytic  solution.  The  analytic  solution  for  the 
slip  was  obtained  from  the  Salvado  and  Minster  solution  described 
Section  2.4.3  of  this  report  by  performing  an  evaluation  of  the  slip 
function  (Equation  (19)  of  Section  2.4.3)  using  exact  Cagniard-de  Hoop 
theory.  The  extraordinarily  good  agreement  confirms  that  the  finite 
difference  method  adequately  treats  the  discontinuity  in  the 
displacement  field  across  the  fault. 
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Joint  Plane 


Geometry  for  two-dimensional  model  of  explosion-induced  block 
motion.  The  complete  radiation  field  consists  of  the  direct 
wave  from  the  explosion  plus  a  scattered  wave  from  slip  on  the 
joint  plane  as  shown. 
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NORMALIZED  SLIP  sh  /<!> 


Analytic 


Finite  Difference 


DIMENSIONLESS  TIME  t8/h 


2345*789 


r  «  1 .Oh 


23456789 


r  »  2.5h 


a  ?  «  ?  f  T  ?  * 


r  *  5. Oh 


Figure  3.  Comparison  of  analytical  and  finite  difference  solutions  for 

the  slip  time  history  on  a  joint  governed  by  the  linear  boundary 
condition,  Equation  (1),  with  $  *  0. 
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2.3.3  Nonlinear  Simulation 


,  A  second  simulation  was  performed  using  the  same  finite 

difference  method,  but  with  the  nonlinear  boundary  condition 
described  in  Section  2.2.  Figure  4a  shows  the  geometry  of  the 
nonlinear  simulation,  as  seen  in  a  horizontal  plane.  The  explosion 
f  Is  represented  by  a  point  source,  with  a  reduced  displacement 

potential  appropriate  for  tuff.  The  P  and  S  wavespeeds  and  density 
are  3140  m/sec,  1580  m/sec,  and  1950  kg/m^,  respectively.  Two 
faults  are  Included,  at  horizontal  distances  of  100  m  and  200  m  from 
t  the  source,  respectively.  Both  faults  dip  away  from  the  source  at 

60*  (so  the  points  nearest  to  the  source  are  actually  at  87  and  173 
meters,  resp^t-1-^ .  The  lithostatlc  pressure  on  the  fault  is 
taken  to  be  :  ,4  &* wrich  corresponds  to  the  overburden  pressure  at 
I  a  depth  if  *£?  m.  Of  course,  since  the  calculation  Is 

two-dimen^i&nau*  depth  variation  of  pressure  can  be  included. 

The  Mfstsi,  yield,  and  elastic  properties  are  approximately 
those  tf  Canute  Mesa  event  MIGHTY  EPIC.  As  shown  in  Figure  4b, 

»  (from  Kipp  and  Kennedy,  1978)  Fault  I  In  our  simulation  corresponds 

to  one  of  the  mapped  faults  near  this  shot,  on  which 
explosion-driven  slip  was  measured  (points  A  and  B  in  Figure  4b). 
This  affords  us  the  opportunity  to  compare  the  numerical  results  to 
»  the  MIGHTY  EPIC  block  motion  observations.  Our  intent,  of  course, 

is  not  to  model  the  details  of  this  particular  event,  but  rather  to 
test  the  plausibility  of  the  simple  dry  friction  boundary  condition 
(Figure  1)  as  a  model  for  explosion-induced  block  motion. 

Figure  5  shows  the  static  slip  computed  on  Faults  I  and  II  in 
the  nonlinear  simulation.  Also  shown  are  the  observed  amplitudes  of 
slip  induced  by  MIGHTY  EPIC  at  sites  A  and  B,  as  given  by  Kipp  and 
Kennedy  (1978).  In  this  figure,  r  represents  distance  from  the 
symmetry  axis  (l.e.,  from  the  point  nearest  the  shot  on  the  plane  of 
the  fault).  At  Site  B,  the  calculated  and  observed  amplitudes  are 
in  very  close  agreement;  at  Site  A,  the  calculated  slip  is  lower 
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Figure  4.  (a)  Geometry  for  the  nonlinear  simulation  of  block  motion. 

(b)  Faults  and  observations  of  explosion-induced  slip  near 
the  event  MIGHTY  EPIC  (after  Kipp  and  Kennedy,  1978).  Solid 
lines  are  mapped  faults;  dashed  lines  show  the  fault  positions 
in  (a)  superimposed  on  the  map  for  comparison.  Solid  triangles 
are  locations  where  slip  was  observed,  and  open  triangles  are 
observation  points  where  no  slip  was  detected.  The  working 
point  and  observations  are  at  approximately  400  m  depth. 
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[iTil 


Simulation 
Fault  II 


Figure  5.  Static  offset  as  a  function  of  radial  distance  for  the  two  faults 
in  the  nonlinear  block  motion  simulation.  Sites  A  and  B  refer  to 
the  MIGHTY  EPIC  observation  points  shown  in  Figure  4b. 
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than  the  observation  by  about  a  factor  of  2.  In  the  simulation, 
significant  slip  continued  on  Fault  I  out  to  a  radial  distance  r  of 
about  300  m.  Calculated  slip  on  Fault  II  Is  everywhere  negligible. 
The  latter  result  Is  In  accordance  with  the  fact  that  no  block 
motion  was  observed  on  faults  located  more  than  about  150  m  from  the 
MIGHTY  EPIC  working  point,  as  Indicated  by  the  open  triangles  In 
Figure  4b. 

Figures  6  and  7  show  slip  and  shear  stress  time  histories, 
respectively  for  Fault  I.  Note  that  no  reversal  of  slip  direction 
occurs.  Note  also  that  beyond  a  range  of  about  150  m  slip  appears 
to  be  mainly  driven  by  dynamic  stress  concentrations  generated  by 
slip  near  the  center  of  the  fault. 

The  general  agreement  shown  in  Figure  5  between  the  simulation 
and  the  data  strongly  suggests  that  our  friction  model  is  an 
adequate  representation  of  the  block  motion  phenomenon.  Since  the 
model  appears  to  be  consistent  with  both  laboratory  and  field  data, 
our  next  step  will  be  to  apply  the  model  to  simulate  block  motion  In 
three  dimensions. 

2.4  AN  EQUIVALENT  SOURCE  ACCOUNTING  FOR  BLOCK  MOTION 

Block  motion  In  the  vicinity  of  a  buried  nuclear  explosion  has 
been  proposed  by  various  authors  as  the  potential  source  of 
anomalous  radiation  often  observed  both  in  the  near  field  and  in  the 
far  field  from  such  sources  (e.g.,  Toksoz  and  Kehrer,  1972;  Masse, 
1981). 

The  fundamental  difference  between  this  mechanism  and  possible 
alternatives  such  as  prestress  relaxation  or  earthquake  triggering 
Is  that  it  Is  a  purely  passive  mechanism.  In  other  words,  there  is 
an  actual  transfer  of  energy  from  the  normal  monopolar  component  of 
the  field  to  the  anomalous  component.  In  contrast,  tectonic 
mechanisms  must  add  energy  to  the  total  radiation  field  through 


17 


srsrews,  science  aho  sorrwAme 


ppage  for  the  nonlinear  simulation. 


systems.  science  and  softwa 


release  of  strain  energy  stored  in  the  medium  prior  to  the 
detonation.  By  using  a  loosely  bonded  Interface  as  a  model  for 
simulating  block  motion,  Salvado  and  Minster  (1980)  have  shown  that 
significant  conversion  of  P  radiation  to  S  radiation  may  occur  if 
the  bonding  across  the  interface  is  weak.  However,  their  solution 
is  based  on  a  first  motion  approximation  valid  only  at  high 
frequency  and  in  the  far  field. 

In  this  report,  we  consider  the  general  problem  of 
axi symmetric  slip  along  a  planar  interface  and  calculate  the 
associated  elastodynamic  radiation.  This  model  is  general  enough  to 
be  appropriate  for  modeling  block  motion  along  a  planar  interface 
near  a  buried  explosion  in  an  otherwise  uniform,  isotropic  infinite 
space.  If  the  slip  function  can  be  determined  as  a  function  of 
position  along  the  interface  for  an  arbitrary  boundary  condition, 
our  solution  is  sufficient  to  calculate  completely  the  anomalous 
contributions  to  the  radiation  field.  In  a  second  section,  we  show 
how  the  general  solution  to  the  problem  attacked  by  Salvado  and 
Minster  (1980)  can  be  obtained  as  a  particular  case,  and  how  a 
rather  simple  expression  can  be  derived  for  the  equivalent  source. 

One  of  the  main  conclusions  of  this  study  is  that  axi symmetric 
block  motion  models  of  this  type  become  poor  r^iators  of  long 
period  energy  under  very  general,  physically  r# a?  Stic  conditions. 
Thus,  in  order  to  explain  observations  of  Rayleigh  wave  phase 
reversals  for  explosive  sources  (Rygg,  1979),  or  observations  of 
strong  Love  wave  radiation  from  such  sources  (e.g.,  Aki  and  Tsai, 
1972),  we  must  either  turn  to  a  different  mechanism,  such  as 
tectonic  release,  or  at  the  very  least,  assume  that  block  motion 
does  not  possess  cylindrical  symmetry.  Departures  from  the 
axlsymmetrlc  geometry  cannot  be  handled  easily  analytically  and 
require  a  numerical  approach. 
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2.4.1  Displacements  from  a  Dislocation  Surface  with  Cylindrical 
Symmetry 

The  displacement  u  caused  by  slip  on  a  surface  in  an  elastic 
medium  is  given  by  the  Volterra  relation: 


“«<&•“>  -  ( 

AUi(x0,«)Cijk£AyGj 

IXjS 

(x,Xo,»)nj(x0)d  Aq 

(2) 

so 

where  a:  is  the  Green's  tensor  for 

the  medium,  Sg  is 

the 

surface  on  which 

slip  occurs,  n  is  the 

normal  to  the  surface. 

x° 

*** 

is  the  coordinate  of  an  integration  point  on  the  surface  and  x  is 

the  observation  point.  For  an  isotropic  medium,  the  elastic  tensor 

C  has  components 
« 
a 

Cijk4  "  X5ij  sk£  +  u^U  4jk  +  4ik  4ji*  * 

For  a  shear  dislocation,  the  slip  au  and  normal  n  vectors  are 
orthogonal.  We  separate  the  slip  into  a  frequency  dependent  scalar 
function  and  a  unit  vector,  au(x°,«)  -  s(x°,<i»)  s  (x°). 

»v  m  cv  ^ 

Equation  (2)  then  simplifies  to 


um(£»w)  *  M  f  +  *  j  "l>  -II  Gjfe0’w)  ^0  * 

J  *xi 
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(3) 


If  we  now  restrict  attention  to  far-field  body  waves.  Equation 
(3)  simplifies  further.  The  far-field  Green's  tensors  may  be 
written  for  P-waves 


G. 

a  J 


A  A  4  1.  r  A 

«  r* .  r  ^  a  ik  xn.r 
m  j  m  e  _  ar~Q 


and  for  S-waves 


A  A 
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fl  j 
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(4a) 


(4b) 
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where  r  Is  a  unit  vector  in  the  direction  of  the  observation  point 
and  ky  •  where  v  -  <*  or  a.  The  coordinates  are  shown  In 
Figure  8.  The  dislocation  Is  located  In  the  z  «  0  plane.  The 
source  position  vector  ^  may  be  written  jtg  •  R  eR  where  R  Is 

the  usual  cylindrical  coordinate  and  eR  is  the  radial  unit 

vector.  We  define  the  cylindrical  coordinate  6  to  be  zero  in  the 

direction  of  the  projection  of  r  on  the  z  »  0  plane.  Note  that  the 
derivative  in  Equation  (3)  simply  brings  down  a  factor  of  ikyr. 

in  the  farfleldi  so  using  1«s  ■  s  Equation  (3)  becomes 


or 


,“.<£•*>  •  » /  5<JSo"“><Vj  *  vG"  “c 


u  *  «-<k»r 


Ik  ju.r 

)  e  ^  dAr 


where 


,A  A  .  A  A  ,  A  A  A 

R“  »  (s.n.  +  s.nj  r,  r.  r 
a  ij  j  1  i  j  m 


(5) 


(6) 


(7) 


.  •  <s1nj  *  sjV  r1  <‘4»  -  rjrm> 


,A  A 


A  A 


A  A 


(8) 


If  we  now  impose  cylindrical  symmetry  on  the  slip  function  so  that 
s  ■  s(R,<i*)  and  s  *  eR,  we  find 

A  A 

s.r  »  sin  e  cos  6 

A  A 

n.r  »  cos  e 
Xg.r  -  R  sin  9  cos  t 

After  some  algebra,  we  find 
R  -  r  sin  2e  cos  6 

Rfl  -  9  cos  29  cos  i 


(An  additional  term  in  R0  as  defined  In  Equation  (8)  Integrates  to 
zero  by  symmetry).  The  displacements  become 
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The  angular  integrals  may  be  performed  using  the  integral 
representation  of  the  iessel  function  (Morse  and  Feshbach,  1953) 

/2ir 

eix  cos*  cos  mi  dt  .  2tim  JJx)  .  (10) 

The  displacements  may  therefore  be  written: 


2.4.2  Generation  of  Long  Period  Body  Waves  and  Surface  Waves 

In  the  long  period  limit  (k  Rsine  «  1)  the  Bessel  function 
becomes 


J,(k  Rsine)  »  1/2  k  Rsine 

•L  3  3 

The  far-field  body  waves  may  then  be  written 
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An  Important  conclusion  may  be  drawn  from  these  equations.  If 
slip  occurs  for  a  finite  time  and  over  a  finite  area,  then  the 
Integrals  in  Equation  (13)  will  approach  a  constant  in  the  long 
period  limit.  The  long  period  amplitudes  of  the  body  waves  are 
therefore  proportional  to  frequency  in  contrast  with  body  waves  from 
an  explosion  or  double  couple  which  approach  a  constant  in  the  long 
period  limit.  This  means  that  the  scattered  radiation  from  an 
explosion  near  a  crack  (excluding  prestress  and  other  nonsymmetrlcal 
effects)  will  not  produce  long  period  Love  and  Rayleigh  waves. 
Axially  synmetrlc  block  motion  of  the  kind  considered  here  is 
therefore  not  likely  to  be  a  viable  explanation  for  anomalous 
surface  waves  sometimes  observed  from  underground  nuclear  explosion. 

The  radiation  patterns  may  be  expressed  in  terms  of  Legendre 
functions;  so  the  displacements  may  be  written  in  terms  of  vector 
spherical  harmonics  of  order  zero.  The  radiation  pattern  factors 
are: 


sin  2e  sin  e  »  £  [P^cos  e)  -  Pj(c os  «)]  , 
d  3  4 

cos  2e  sin  e  «  [|-  P^cos  e)  -  P3(cos  e)] 


(H) 


The  radiation  patterns  contain  only  Legendre  functions  of  degree  1 
and  3  corresponding  to  a  dipole  term  plus  an  octupole  term. 

In  terms  of  vector  spherical  harmonics  (e.g.,  Stevens  1980), 
the  displacement  may  be  written 


U(x,«)  -  S  b£0(«)  U£0  (kBr)  +  Yg0(«)  L£0  (kar) 
£  .  1,3 
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and  where 

/op 

dR  R2  S(«,R)  . 

u 

Equation  (15)  is  valid  in  the  near-field,  as  well  as  in  the 
far-field,  up  to  the  smallest  sphere  enclosing  the  region  of  slip  at 
all  periods  long  enough  that  Equation  (12)  remains  valid. 
Near-field  waveforms  and  surface  .waves  may  also  be  computed  using 
the  coefficients  in  Equation  (15). 

The  radiation  field  may  also  be  expressed  in  terms  of  the 
scalar  Cartesian  potentials  defined  by 


X1  -  \  (vxu)i  i  -  1,2,3 

a 

X  «Y-u 


These  potentials  are  harmonic  functions  and  may  be  expanded  in 
scalar  spherical  harmonics: 

xV.e.P.u,)  «S  2  (A1  (M)cos  +  B.1  («)sin  Pem(cos  e)  hfl2)(kvr) 

which  is  a  useful  form  for  synthesizing  body  waves  and  surface  waves 
in  layered  media  (e.g.,  Bache  and  Harkrider,  1976). 

The  nonvanishing  multipole  coefficients  are: 
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Given  the  slip  function  as  a  function  of  position  along  the 
interface  across  which  the  dislocation  occurs,  we  can  compute  I,  and 
hence,  the  complete  radiation  field.  In  the  next  section,  we  apply 
this  formalism  to  the  case  of  block  motion  triggered  by  the 
detonation  of  an  explosion  near  a  linearly  viscous  interface. 

2.4.3  Seismic  Radiation  from  Explosion  Near  a  Viscous  Interface 

The  problem  of  scattered  radiation  from  an  explosion  near  a 
viscous  plane  interface  was  treated  by  Salvado  and  Minster  (1980) 
using  a  Cagniard-de  Hoop  first  motion  approximation.  Using  the  slip 
function  derived  by  Salvado  and  Minster  and  the  equations  derived  in 
the  last  section,  we  can  obtain  exact  solutions  for  the  scattered 
field  in  the  frequency  domain. 

The  potential  for  an  explosion  in  an  infinite  space  may  be 
written 


A<4) 

*30 

A<2> 

A11 

A<2> 

A31 
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where  y(u)  is  the  reduced  displacement  potential  and  displacement  is 
given  by 


ur(r,«) 
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In  cylindrical  coordinates  this  may  be  written 


«(R,z,«)  -  -  'F(u)  / 

A  •/ 

*~0(kR) 

“vai*“h|  k 

e  *  —  dk 

V 

a 

(18) 

Q 

where 

2  2 
v  -  (r  -  k2) 
a  '  o' 

k  >  k 

a 

2  2 

-  i(k2  -  k2). 

k.<  k 

a 

and  the  explosion  is  located  at  R 

«  0,  z  -  h  (see  Figure  9) . 

The  slip  function 

s  on 

the  interface  is  defined 

by 

u(z  -  0+)  -  u(z  »  0”). 

Using 

the  notation  of  Salvado 

and 

Minster  and  their  Equations  (2.7)  -  (2.8)  (correcting  a  sign  error 
in  (2.8a)),  the  slip  velocity  s  is  found  to  be 


/°°  aUn 

dk  2[vfl  Sl  (k)  + 


(k)] 


(19) 


where 


jm  OQ 

-  V  («)  4up(1  -  4)  J  dk 

<r 


n  if  —V  h 

ICO^kR)  \>B  u  e  a 


W 


2upv  iw^ 


9  9  «•  9  9  ‘•WPV 

D(k)  ■  [(2wk  -  ou  )  -  4u  v  v  k2]  (1  -  $)  - 

P 


Here  4>  is  a  parameter  which  represents  the  amount  of  viscous 
damping.  4>  -  1  implies  perfect  bonding  while  $  »  0  implies  an 
interface  with  no  bonding.  v„  is  defined  similarly  to  v  but 
with  a  replaced  by  the  shear  velocity  s. 


Substituting  this  expression  into  the  surface 
(Equation  (11)),  the  Integrals  have  the  form 


./ 


dR  R  Jx(kB  R  sin  e)  f  dk  k  f(k)  Jj(kR) 

<T 


integrals 
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Plane  Viscous  Interface 


Figure  9 


Coordinates  used  for  explosion  near  a  viscous  interface. 
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Interchanging  the  order  of  integration  and  using  the  identity 

f  ~  Jm0cR)  Om(k*R)  dR  -  --Hi  , 

ty 

the  integrals  vanish  leaving 
I  -  f(k) 


where  k  ■  ^  sin  e  where  v  «  a  for  P-waves  and  8  for  shear  waves. 
The  displacements  simplify  to  the  following  form 


uy(r,«»“) 


1 


y«> 


-«n  h 
•  ® 


(20) 


The  last  factor  in  this  expression  is  the  spectrum  of  the 
direct  P-wave  from  the  explosion  (except  for  a  phase  factor  due  to 
propagation).  ky  «  u/v;  ey  *  r  for  P-waves  and  e  for  S-waves. 
We  define  the  following  quantities. 

For  P-waves: 

k  «  sin  e/a 

i  2  if  i 

na  -  7  U  -  sin  e)  -  -  {cos  e| 

i  a2  2  1/2 

aa  *  ?  (1  ”  T  sin  9) 

a 

Ay(e)  m  sin  2e 


For  S-waves: 

k  «  sin  e/6 


i  2  1/2  i  t  I 

nB  «  j  (1  -  sin^  e)  -  j  [cos  e| 

i  a2  2  1/2 

n  ■  -  (1  -  •=T  sin  e)  if  sin  e  <  8/a 

a  a  ac  — 
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if  sin  e  >  a/a 
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Then  Ty(e)  is  a  frequency  independent  function  defined  by 
[2ykoin  aa2/v3](l  -  *  ) 

V>  -  - hsm - V> 


(21) 


where 


D(k)  «  ({2yk2  -  p)  -4y2na  ng  k2)(l  -4  )  -  2yping  4/b 


The  angular  function  Ty(e)  is  real  for  all  angles  for  the 
P-wave,  and  for  angles  less  than  the  critical  angle 
(e  ■  sin-*  s/a)  for  the  S-wave.  It  is  complex  for  angles  in 
the  range  e  <  e  <  *  -  e  .  For  P-waves  and  shear  waves  outside 
the  critical  range,  the  waveforms  are  identical  to  the  explosion 
waveform  except  for  an  angle  dependent  amplitude  factor.  Within  the 
critical  range,  the  waveforms  are  phase  shifted  because  of  the 
complex  value  of  Tfl(e),  and  also  damped  at  high  frequency  by  the 
exponential  function  e-^a11  . 

The  radiation  pattern  functions  Ty(e)  are  shown  for  several 
values  of  the  damping  parameter  $  in  Figures  10  (P-waves)  and  11 
(S-waves).  These  radiation  patterns  do  not  agree  with  the  shear 
wave  radiation  patterns  of  Salvado  and  Minster  (their  Figure  4)  for 
angles  beyond  the  critical  angle.  Their  equations  (Table  1)  do, 
however,  agree  with  our  equations  for  the  radiation  pattern;  so 
apparently  their  figure  is  in  error.  The  radiation  patterns  found 
here  for  4-0  are  also  identical  to  radiation  patterns  for  the 
waves  from  a  strike  slip  double  couple  reflected  from  a  free  surface 
found  by  Burridge,  Lapwood  and  Knopoff  (1964). 
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Figure  10.  Far-field  P-wave  radiation  pattern  amplitudes  jT (0 ) |  from 
an  explosion  near  a  viscous  Interface  for  dampingaparameters 
of  0.0,  0.1,  O.S,  and  0.9.  The  maximum  amplitudes  on  the  figure 
are  0.47,  0.39,  0.12,  and  0.14  respectively.  Ta(e)  is  real 
at  all  angles,  positive  In  the  upper  half-space,  negative  in 
the  lower  half-space. 
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Figure  11  • 


Far-fleld  shear  wave  radiation  pattern  amplitudes  |Tg (8) |  from 
an  explosion  near  a  viscous  interface.  Figures  are  for  damping 
parameters  of  0.0,  0.1,  0.5  and  0.9.  The  maximum  amplitudes  on 
the  figures  are  5.2,  1.8,  0.54,  and  0.14  respectively.  Beyond 
the  critical  angle  Te(Q )  is  complex  and  the  amplitudes  are  ex¬ 
ponentially  damped  with  frequency.  Tg(9)  is  real  and  positive 
outside  the  critical  range. 
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In  these  examples  we  have  used  a  -  6.0,  s  -  3.527,  o  «  2.0;  so 
the  critical  angle  accurs  at  e£  -  36  degrees.  The  sharp  spikes 
for  small  values  of  $  occur  at  the  critical  angles  and  have  a 
maximum  value  of  a/s  tan  2e  as  5.  For  the  unbonded  case, 

V 

therefore,  the  scattered  shear  wave  will  exceed  the  direct  P-wave  in 
amplitude. 

It  is  interesting  to  note  that  in  the  long  period  limit,  the 
displacements  approach  a  non-zero  constant.  Earlier,  we  showed  that 
for  any  slip  which  occurs  over  a  finite  region  and  for  a  finite 
time,  the  spectra  should  be  proportional  to  frequency  at  long 
periods.  The  reason  for  this  anomalous  behavior  is  the  linearity  of 
the  boundary  condition.  This  allows  the  slip  on  the  boundary  to 
propagate  to  infinity  and  changes  the  long  period  content  of  the 
waveform.  The  slip  function  given  in  Equation  (19)  for  the  case 
<&  ■  0  is,  in  fact,  identical  to  the  surface  radial  velocity  caused 
by  an  explosion  in  a  uniform  half-space,  a  rather  surprising 
conclusion,  considering  the  difference  in  continuity  conditions 
across  the  boundary. 


2.4.4  Analysis  of  Scattered  Wave  from  Finite  Difference  Calculation 

Equation  (13)  gives  the  scattered  radiation  at  long  periods 
due  to  block  motion  with  cylindrical  synroetry.  In  the  long  period 
limit,  the  slip  velocity  is  equal  to  the  static  limit  of  the  slip 
s(t  «*)  =  s*.  We  can  define  the  constant 


4f< 

2 a  J 


dR  IT  sjR) 


Then,  at  sufficiently  long  periods,  the  scattered  P-wave  plus  the 
direct  P-wave  is  given  by 


u 


-Ik  r 

'  *  1f  I«s1n2esine] 
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where  is  the  static  limit  of  the  explosion  reduced  velocity 
potential  and  f  is  frequency. 

Uo  Mas  found  by  performing  a  numerical  integration  of  the 

static  limit  of  the  slip  function  for  the  two  nonlinear  finite 
difference  models.  In  the  first  model 

I*  -  .032  ^ 

while  In  the  second  model 
I*  -  *<#3 

Clearly,  the  scattered  wave  Is  very  small  compared  to  the 

explosion.  The  approximation  used  Is  valid  to  about  one  hertz.  At 
one  hertz,  the  scattered  wave  is  less  than  five  percent  of  the 

direct  P-wave.  At  20  seconds  it  is  reduced  to  less  than  0.2  percent 

of  the  direct  wave.  Clearly,  this  type  of  scattering  cannot  cause 
Rayleigh  wave  reversals. 


2.5  SUMMARY  AND  CONCLUSIONS 

The  principal  conclusions  of  our  study  of  block  motions  near  a 
burled  explosive  source  can  be  summarized  as  follows. 

1.  Qualitative  effects  of  block  motions  in  axisymmetric 
geometry  on  the  outgoing  radiation  field  are  usefully 
described  by  an  analytical  model  if  the  boundary  condition 
across  the  joint  is  linear.  Comparison  of  analytical 
results  with  finite  difference  axisymmetric  calculations 
showed  that  our  numerical  treatment  of  Imperfectly  bonded 
joints  in  the  rock  mass  leads  to  correct  and  accurate 
results. 

2.  Quantitative  agreement  between  numerical  simulations  and 
field  observations  is  achieved  for  a  rate  independent  dry 
friction  boundary  condition  with  a  kinetic  friction 
coefficient  of  0.77. 
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3.  In  axisynmetric  geometry  an  equivalent  seismic  source 
which  account  for  block  motion  contains  both  dipole  and 
octupole  terms.  Evaluation  of  these  contributions  to  the 
radiation  field  for  both  linear  and  nonlinear  boundary 
condition  shows  that  this  equivalent  source  is  a  very 
poor  radiator  at  long  periods,  and  cannot  lead  to 
Rayleigh  wave  phase  reversals.  It  could  affect  body  wave 
radiation  significantly,  however. 

8ased  on  these  results,  it  appears  that  three-dimensional 
effects  which  lead  to  departures  from  axisymmetry  should  be  studied 
next.  The  most  likely  situation  In  which  our  conclusions  to  date 
might  be  invalid  is  one  where  block  motion  occurs  along  a  45* 
slip-slip  planar  joint  with  strong  asymmetry.  We  are  currently 
investigating  this  case  through  a  fully  three-dimensional  numerical 
calculation. 
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III.  EARTHQUAKE  COUPLING  TO  ACOUSTIC  WAVES 


3.1  INTRODUCTION 

In  this  section  we  summarize  our  progress  in  a  study  to 
simulate  the  acoustic  waves  generated  by  small  earthquakes  and 
explosions.  S3's  part  of  this  effort  is  to  provide  seismological 
input  to  the  Mission  Research  Corporation  (MRC)  who  will  perform  the 
acoustic  wave  simulations. 

There  are  four  tasks  to  be  performed  under  this  contract.  The 
first  task  is  a  validation  study  of  previously  performed  sir»u1ations 
for  the  1971  San  Fernando,  California,  earthquake.  Assumptions  used 
in  the  previous  calculations  will  be  checked  by  comparison  of  the 
previous  results  to  more  rigorous  simulations  for  selected  cases. 
The  second  task  consists  of  a  review  of  the  seismic  literature  and 
near-field  ground  motion  data  base  to  determine  the  range  of  source 
characteristics  possible  for  earthquakes  of  approximately  magnitude 
four  and  to  choose  a  representative  event  for  modeling.  The  third 
task  Is  to  obtain  a  source  model  for  the  chosen  event  which  can  be 
used  to  calculate  the  near-field  vertical  ground  accelerations. 
These  will  be  used  by  MRC  for  simulation  of  the  acoustic  response. 
The  final  task  is  to  provide  our  earthquake  modeling  results  to  MRC, 
together  with  other  seismological  input  needed  to  perform 
simulations  of  the  acoustic  response  to  small  magnitude  earthquakes 
and  explosions. 

In  the  sections  to  follow  we  will  sunmarize  our  progress  in 
Tasks  1  through  3. 

3.2  VALIDATION  OF  PREVIOUS  ACOUSTIC  WAVE  SIMULATIONS 

Under  a  previous  contract  we  performed  a  simulation  of  the 
near-field  ground  motions  from  the  1971  San  Fernando,  California, 
earthquake  and  provided  the  results  to  MRC  for  estimation  of  the 
acoustic  response  to  that  event.  The  results  of  this  study  have 
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been  reported  by  Bache,  et  al_.  (1981).  For  each  time  point  in  the 
acoustic  signature,  a  surface  integral  of  surface  ground 
acceleration  must  be  calculated.  Computational  considerations 
demand  a  simple  and  inexpensive  algorithm  for  generating  the 
near-field  ground  motions.  To  satisfy  this  need,  we  have  developed 
rather  simple  models  of  the  seismic  energy  radiation  from  the 
earthquake  and  of  the  interaction  of  this  radiation  field  with  the 
free  surface.  Several  approximations  were  necessary  and  the  purpose 
of  this  task  Is  to  examine  the  validity  of  these  approximations. 

The  free  surface  response  used  is  a  high  frequency 
approximation  first  proposed  by  Knopoff,  et  al_.  (1957).  It  is 
reasonably  accurate  near  the  geometric  arrivals  but  contains  no 
Rayleigh  waves  or  S-wave  to  P-wave  converted  phases.  When  using 
this  approximate  free-surface  response,  the  acoustic  signatures 
generated  appear  to  be  consistent  with  those  estimated  using 
geometrical  optics  (W.  Wortman,  personal  communications).  In  other 
words,  the  acoustic  signals  in  the  linear  portion  of  the  atmosphere 
have  a  character  very  much  like  the  vertical  ground  velocity  at  a 
point  on  the  free  surface  along  the  geometrical  ray  path  from  source 
to  observer. 

Calculations  have  been  performed  to  determine  whether  these 
results  still  hold  if  a  rigorous  free  surface  response  is  employed. 
Green's  functions  for  the  acoustic  pressure  due  to  a  buried 
dislocation  were  computed  using  Cagnlard-de  Hoop  techniques.  These 
Green's  functions  Include  all  diffraction  effects  and  converted 
phases.  Calculations  performed  to  date  suggest  that  geometrical 
optics  performs  quite  well  at  least  out  to  ranges  of  a  few  source 
depths.  Further  calculations  are  needed  to  determine  to  what  ranges 
our  free-surface  approximation  remains  adequate. 

Mr.  William  Wortman  of  MRC  has  provided  us  with  a  sample 
calculation  of  an  atmospheric  pressure  signature  derived  from  our 
model  of  the  San  Fernando  earthquake.  We  will  perform  an 
Independent  calculation  of  that  response  using  our  source  model  and 
Green's  function  generated  by  Cagnlard-de  Hoop  methods. 
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3.3  REVIEW  OF  SEISMIC  LITERATURE  AND  NEAR-FIELD  DATA  BASE 

We  have  surveyed  the  seismic  literature  to  determine  the  kinds 
of  source  parameters  to  be  expected  for  earthquakes  of  the  magnitude 
range  of  Interest.  Until  recently,  nearly  all  of  the  Information 
about  such  events  has  come  from  spectral  measurements  and  corner 
frequency  analysis.  Most  notable  are  the  studies  of  Thatcher  and 
Hanks  (1973)  who  examined  the  recordings  of  a  number  of  southern 
California  events,  and  Tucker  and  Brune  (1973),  who  studied  the 
aftershocks  of  the  1971  San  Fernando  earthquake.  These  studies  are 
useful  for  obtaining  an  estimate  of  the  range  of  moments  and  stress 
drops  one  might  expect  for  events  within  the  magnitude  range  of 
Interest.  Stress  drop  may  be  a  very  Important  parameter  controlling 
near-field  ground  motion  characteristics,  particularly  when  the 
corner-frequency  is  very  high,  since  local  magnitudes  or  body  wave 
magnitudes  are  not  sensitive  to  the  additional  energy  present  at 
high  frequencies.  It  Is  clear  from  these  studies  that  the  range  of 
possible  stress  drops  may  span  two  orders  of  magnitude  with  an 
average  of  about  50  to  100  bars. 

Detailed  modeling  of  small  magnitude  earthquakes  have  recently 
been  undertaken  at  the  United  States  Geological  Survey  (USGS),  Menlo 
Park,  California.  Three  aftershock  sequences,  the  1975  Oroville 
earthquakes,  the  1979  Mammoth  Lakes  earthquakes,  and  1979  Livermore 
earthquakes,  have  provided  enough  near-field  data  to  make  detailed 
modeling  possible.  Results  of  these  modeling  studies  will  soon  be 
available  in  USGS  open  file  reports. 

For  this  study  the  most  useful  data  set  appears  to  be  the 
Oroville  aftershock  strong  motion  recordings  (Fletcher,  et.  a!., 
1980).  There  are  several  events  In  the  magnitude  range  of  Interest 
which  have  triggered  more  than  five  accelerograph  stations.  All 
data  consist  of  three-component  recordings  with  band-width  of  0.5  to 
20.0  Hz.  Locations  and  fault  plane  solutions  are  available  for  all 
events  and  estimates  of  the  basic  source  parameters  have  also  been 
obtained  for  most.  Stress  drops  tend  to  be  of  the  order  of  100-200 
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bars  (J.  Boatwright,  personal  communication) .  We  have  chosen  one  of 
these  events  for  a  detailed  modeling  study  and  acquired  the  data 
from  the  U.S.  Geological  Survey. 

3.4  MODELING  OF  OROVILLE  AFTERSHOCK 

The  event  we  have  chosen  for  detailed  source  modeling  is  event 
0247  (August  3,  1975)  of  the  Oroville  aftershock  sequence.  Its 

source  mechanism  shows  normal  faulting  at  a  depth  of  7.4  km  and  its 

local  magnitude  (M^)  is  4.1.  Eight  strong-motion  accelerograph.s 

triggered  on  the  event  and  seven  of  the  recordings  are  of  a  quality 
good  enough  for  detailed  modeling.  Figure  12  shows  the  epicenter 
and  the  location  of  the  accelerograph  stations. 

One  of  the  processed  accelerograms  is  shown  in  Figure  13. 

Because  these  data  were  recorded  by  triggered  instruments,  the  first 
arrivals  are  lost,  but  the  shear  wave  arrivals  are  quite  clean  and 
have  rather  simple  waveforms.  We  will  use  ground  velocity  for  our 
modeling  study,  since  It  is  velocity  which  most  closely  reflects  the 
characteristics  of  the  acoustic  signals  we  eventually  wish  to 
simulate.  A  processor  for  computing  ground  velocity  from  the 
accelerograms  has  been  kindly  provided  with  the  data  by 
J.  Boatwright  of  the  USGS  A  sample  of  computed  velocity  is  shown  in 
Figure  13.  Estimates  of  source  parameters  have  also  been  provided 
by  the  USGS 

Our  approach  is  to  model  the  ground  velocity  with  the  methods 
described  in  Bache,  et_  (1981).  A  source  representation  similar 
to  that  which  we  use  here  has  already  been  successfully  employed  for 
other  Oroville  events  (Boatwright,  1981).  Once  the  modeling  is 
completed,  the  results  will  be  forwarded  to  MRC. 
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Figure  12.  Epi central  region  of  Oroville  aftershock 
0247,  showing  the  location  of  nearby 
recording  stations  (after  Boatwright, 
1981).  The  depth  of  the  event  is  7.4  km. 
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